First study of the gluon-quark-antiquark static potential in SU(3) Lattice QCD 
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We study the long distance interaction for hybrid hadrons, with a static gluon, a quark and an 
antiquark with lattice QCD techniques. A Wilson loop adequate to the static hybrid three-body 
system is developed and, using a 24 3 x 48 periodic lattice with (3 = 6.2 and a ~ 0.075 fm, two 
different geometries for the gluon-quark segment and the gluon-antiquark segment are investigated. 
When these segments are perpendicular, the static potential is compatible with confinement realized 
with a pair of fundamental strings, one linking the gluon to the quark and another linking the same 
gluon to the antiquark. When the segments are parallel and superposed, the total string tension is 
larger and agrees with the Casimir Scaling measured by Bali. This can be interpreted with a type-II 
superconductor analogy for the confinement in QCD, with repulsion of the fundamental strings and 
with the string tension of the first topological excitation of the string (the adjoint string) larger than 
the double of the fundamental string tension. 



I. INTRODUCTION 

Here we explore the static potential of the hybrid three- 
body system composed of a gluon, a quark and an an- 
tiquark using lattice QCD methods. The Wilson loop 
method was deviced to extract from pure-gauge QCD 
the static potential for constituent quarks and to provide 
detailed information on the confinement in QCD. In what 
concerns gluon interactions, the first lattice studies were 
performed by Michael [l], and Bali extended them to 
other SU(3) representations @. Recently Okiharu and 
colleagues [J, § extended the Wilson loop for tree-quark 
baryons to tetraquarks and to pentaquarks. Our study of 
hybrids continues the lattice QCD mapping of the static 
potentials for exotic hadrons. 

The interest in hybrid three-body gluon-quark- 
antiquark systems is increasing in anticipation to the fu- 
ture experiments BESIII at IHEP in Bcijin, GLUEX at 
JLab and PANDA at GSI in Darmstadt, dedicated to 
study the mass range of the charmonium, with a focus 
in its plausible hybrid excitations. Moreover, several ev- 
idences of a gluon effective mass of 600-1000 MeV from 
the Lattice QCD gluon propagator in Landau gauge, 
@j 0] j from Schwinger-Dyson and Bogoliubov-Valatin so- 
lutions for the gluon propagator in Landau gauge Q, 
from the analogy of confinement in QCD to supercondu- 
tivity [9|, from the lattice QCD breaking of the adjoint 
string [1|, from the lattice QCD gluonic excitations of 
the fundamental string [To| from constituent gluon mod- 
els (ill . [l2l [l3l compatible with the lattice QCD glueball 
spectra [l4ll5l. liH Uj\ , and with the Pomeron trajectory 
for high energy scattering [H, Qj| may be suggesting that 
the static interaction for gluons is relevant. 

Importantly, an open question has been residing in the 
potential for hybrid system, where the gluon is a colour 
octet, and where the quark and antiquark are combined 
to produce a second colour octet. While the constituent 
quark (antiquark) is usually assumed to couple to a fun- 
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FIG. 1: String attraction and fusion, and string repulsion, 
respectively in type I and II superconductors 



damental string, in constituent gluon models the con- 
stituent gluon is usually assumed to couple to an adjoint 
string. Notice that in lattice QCD, using the adjoint 
representation of SU(3), Bali [3j found that the adjoint 
string is compatible with the Casimir scaling, were the 
Casimir invariant Xi ■ Xj produces a factor of 9/4 from 
the qq interaction to the gg interaction. Thus we already 
know that the string tension, or energy per unit lenght, 
of the adjoint string is 1.125 times larger than the sum 
of the string tension of two fundamental strings. How 
can these two pictures, of one adjoint string and of two 
fundamental strings, with different total string tensions, 
match? This question is also related to the superconduc- 
tivity model for confinement, is QCD similar to a Type-I 
or Type-II superconductor? Notice that in type Type- 
II superconductors the flux tubes repel each other while 
in Type-I superconductors they attract each other and 
tend to fuse in excited vortices. This is sketched in Fig. 
[TJ The understanding of the hybrid potential will answer 
these questions. 

In Section II we produce a Wilson Loop adequate 
to study the static hybrid potential. In Section III wc 
present the results of our Monte-Carlo simulation, in a 
24 3 x 48 pure gauge lattice for (3 = 6.2, corresponding to 
a lattice size of (1.74 fm) 3 x (3.48 fm), assuming a string 
tension ^fa = 440 MeV. In Section IV we interpret the 
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FIG. 2: Wilson loop for the gqq potential 



results and conclude. 



II. HYBRID WILSON LOOP 

In principle any Wilson loop with a geometry similar 
to the one in Fig. [2] and describing correctly the quantum 
numbers of the hybrid is adequate, although the signal 
to noise ratio may depend in the choice of the Wilson 
loop. A correct Wilson loop must include a SU(3) octet, 
the gluon, a SU(3) triplet, the quark and a SU(3) anti- 
triplet, the antiquark. It must also include the connection 
between the three links of the gluon, the quark and the 
antiquark. 

In the limit of infinite quark mass, a nonrelativistic po- 
tential V can be derived from the large time behaviour of 
cuclidean time propagators. Typically, one has a meson 
operator O and computes the Green function, 



O(t)O(0)|0) 



exp{-Vt} 



(1) 



for large t. Different types of operators allow the defi- 
nition of different potentials. In the static gluon-quark- 
antiquark interaction, the static gluon can be replaced 
by a static quark-antiquark pair in a colour octet rep- 
resentation. In this way, we can construct the gluon- 
quark-antiquark Wilson loop starting from the mesonic 
operator, 



0(x) 



1 



q{x) X a r x q(x) q{x) X a T 2 q{x 



(2) 



where Ti are spinor matrices. Using the lattice links to 
comply with gauge invariance, the second operator in cq. 
@ can be made nonlocal to separate the quark and the 
antiquark from the gluon, 



0{x) 



q(x) X a T iq (x) 



q(x - rifii)!/^ (x - nfix) ■ ■ ■ (x - /ti) 
A Q T 2 E/ M2 (z)---t/ /l2 (x+(r 2 -l)£2) 
q{x + r 2 jl 2 ) ■ (3) 



The nonrelativisit potential requires the computation of 
the Green functions present in eq. {!]). Assuming that 
all quarks are of different nature, the contraction of the 
quark field operators gives rise to the gluon operator, 

rlrhllit - l,x) ■ ■ -Ul(0,x) X b 



Wo = 
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U 4 (0,x) 



Tr 



{ U^ 2 (t,x) 



Ui{t-l,x) A a | x 
• U tl2 (t,x+ (r 2 - l)/t2) 



U\{t -l,x + r 2 fi2) ■ ■ ■ Ut(0, x + r 2 fe) 
Ul 2 (0,x + (r 2 -l)fi 2 )---Ul 2 (0,x) X" 

U't i (0,ar-Ai)"-^ 1 (0 J a!-r 1 Ai) 
C/4(0, x - nfbi) ■ ■ ■ Ui{t - 1, x - ri/ti) 

U^{t,x-n(n)---U^{t,x- [ix) A a j . (4) 

Gauge invariance of ((J)) can be proven with the help of 
the relation 

\ a \ \ \ 

—J = ^Sirfjk - ^SijSki . (5) 
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How docs our operator relate with the operators used 
so far to investigate the gluon interactions on the lattice? 
The gluonic time-like links used by Michael and collogues 
[2I. [20(1 to study the glue lump are the real 8x8 matrices, 



^Adj a/3 = }_ Tr{UiX a UlX 13 } 



(0) 



built from the usual SU(3) fundamental representation 
links Ui jwhereas in the investigation of Casimir scaling 
by Bali [3J], the author worked directly with adjoint links, 
i.e. with the 8x8 matrix SU(3) representation. If one 
now compares the Wilson loop in cq. ([4]) with cq. ([6]), it 
comes that when t corresponds to a single lattice spacing, 
then the gluonic trace, i. e. the first trace, of eq. (J4j) is 
a "Michael link" . 



III. THE STATIC HYBRID POTENTIAL 

In this paper we consider two possible hybrid geome- 
tries: _L with the quark-gluon segment perpendicular to 
the antiquark-gluon segment; || with the quark-gluon seg- 
ment parallel to the antiquark-gluon segment. We denote 
the potentials, respectively, V±(ri,r 2 ) and V\\(ri,r 2 ), 
where r\ (r 2 ) is the quark-gluon (antiquark-gluon) dis- 
tance in lattice units, defined in eq. (|3]). 

Here we discuss the results of our simulation with 142 
24 3 x 48, /3 = 6.2 pure-gauge Wilson action SU(3) con- 
figurations. The configurations are generated with the 
version 6 of the MILC code [2l|, via a combination of 
Cabbibo-Mariani and overrelaxed updates. In order to 
improve the signal to noise ratio, the links are replaced 
by "fat links" [22| 

UpW - YT6^ (s) + w ^ u ^ u ^ s + ( B + A 
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FIG. 3: Potential in lattice spacing units, for the system qqg, 
respectively for the _L and || geometry 

+U}(a-v)Ur(a-v)U v (a + n-v)) (7) 

followed by a projection into SU(3). We use w = 0.2 and 
iterate this "fat link" 25 times both in the spatial direc- 
tion and in the temporal one. The temporal smearing 
reduces the short-range Coulomb potential but produces 
a clearer signal for the long-range potential, the one we 
are interested in. Furthermore, to improve the quality 
of the signal, we explore the symmetry T\ <-» r 2 when 
computing V±(ri,r 2 ) and Vj^fi,^). 

Using eq. ([I]), the static potentials are extracted from 
the fit of minus the log of the Wilson loop, — log W, for 
large euclidian time t. This fit provides us with the po- 
tential, and we estimate the respective error bar with the 
jackknifc method. 

We are essentially interested in the largest possible dis- 
tances, to compare the different possible string tensions. 
With 24 3 x 48 periodic lattices, the maximum distance we 
reach is 10 a. In this way we avoid possible finite volume 
effects. At 12 a, due to our periodic boundary conditions, 
we already could approach a maximum, deviating from 
the linear behaviour. Notice that we calibrate our lat- 
tice spacing a ~ 0.075 fm as in Bali and Schilling [23[. 
Thus our maximal distance is still comfortably shorter 
than the string breaking distance, larger than 1 fm, and 
comfortably longer than the pertubative distance of say, 
0.3 fm. We also start measuring the potentials at the 
distance of 6 a because we are interested in studying the 
long-distance, non-pcrtubative part of the static poten- 
tials. Our results for the static hybrid potentials V± and 
Vj| are displayed in Fig. [3] 

Again, to get the string tensions cr, we fit the large 



FIG. 4: Potential for the system qqg, in lattice spacing units, 
respectively for the _L geometry as a function of n + ri and 
for the || geometry as a function of n = ri — r 

distance part of the potentials with a linear potential, 
with the same method we used for the temporal fit to ex- 
tract the static potentials. We admitting a singlet quark- 
antiquark singlet string tension of 440 MeV, which cor- 
responds to an inverse lattice spacing of a -1 = 2718 ± 32 
MeV, according to Bali and Schilling (23|. 

To study the onset of two fundamental strings, we plot 
in Fig. 2] the perpendicular geometry potential V±_ as a 
function of the sum of the two distances in lattice spacing 
units, r\ between the quark and the gluon and r 2 between 
the antiquark and the gluon, as in cq. ([!]). Indeed the 
potential is linear in the sum of the distances. We further 
fit the potential to V±(ri,T2) = Co + cr(ri + T2), and get 

sfd= 441 ±6 MeV = (1.00 ± 0.01)^ (8) 

(x 2 /dof = 1.34802) which is consistent with a = a a , 
reinforcing the picture that, at long distances we have 
two fundamental strings one linking the quark and the 
gluon and, the other, linking the antiquark to the gluon. 

To study if the double fundamental string picture can 
also be compatible with the Casimir scaling result found 
by Bali 0, we also consider the case where the quark and 
antiquark are superposed. In this case the static quark 
and antiquark are equivalent to a static gluon, and there- 
fore our potential is equivalent to a static gluon-gluon po- 
tential. This is the case of the parallel geometry potential 
Vy when the two distances t\ between the quark and the 
gluon and r 2 between the antiquark and the gluon, as in 
cq. ((3]), are identical, r\ = r 2 . This is plotted in Fig. 0] 
and indeed we find a linear behaviour. We further fit the 
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static potential of the parallel geometry with r\ = ti = r, 
by Vj| (r, r) = c + or, and we get 

N /a = 681±9MeV= (1.55 ± 0.02)^/^0 (9) 

(x 2 /dof = 0.235172) wich is a little larger than the 
Casimir scaling y/a = fy/oo" ratio. Nevertheless our re- 
sults confirm with the ones of Bali Q since our error 
bars match. The increase or the static hybrid poten- 
tial when the quark and antiquark are superposed can 
be interpreted with a repulsive energy between the two 
fundamental strings. This repulsive energy also exists in 
Type-II superconductors. 

IV. CONCLUSION 

We explore, in 24 3 x 48 periodic lattices with (3 = 6.2 
and a ~ 0.075, two different geometries for the gluon- 
quark segment and the gluon-antiquark segment. When 
these segments are perpendicular, the static potential is 
consistent with confinement realized with a pair of fun- 
damental strings, one linking the gluon to the quark link- 
ing the same gluon to the antiquark. When the segments 
are parallel and superposed, the total string tension is 
larger and is compatible with a repulsive energy between 
the two fundamental strings. Notice that when the two 
segments are parallel and superposed, the total string 
tension is also compatible with the Casimir Scaling mea- 
sured by Bali. 

This can be interpreted with a type-II superconductor 
analogy for the confinement in QCD, with repulsion of 



the fundamental strings and with the string tension of 
the first topological excitation of the string (the adjoint 
string) larger than the double of the fundamental string 
tension. Nevertheless, because the energy of two funda- 
mental strings plus the repulsive energy measured here 
is quite similar to the energy of the adjoint string mea- 
sur ed by Bali 0, this shows that the pure gauge QCD 
is similar to a Type-II superconductor quite close to the 
phase transition to a Type-I superconductor. 

Our results are important for constituent models for 
hybrids and glueballs. In the three-body hybrid, with one 
quark, one antiquark and one gluon, our results suggest 
that the best potential model has only two fundamen- 
tal strings, plus a repulsion acting only when the two 
fundamental strings are close. In the two body gluon- 
gluon glucball, our results suggest that the string ten- 
sion is similar to the one of the Casimir Scaling model, 
with a factor of the order of | when compared with the 
quark-antiquark potential. We can also extrapolate our 
result for three-body glueballs, relevant for the odderon 
problem. With three gluons, a triangle formed by three 
fundamental strings costs less energy than three adjoint 
strings with a starfish-like geometry. Thus we antici- 
pate that the thrcc-gluon potential is similar to a sum 
of three mesonic quark-antiquark interactions, plus a re- 
pulsion acting only when there is superposition of the 
fundamental strings. 
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